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We present a complete treatment of the diffusion processes for supersymmetric 
electrowealc baryogenesis tliat characterizes transport dynamics ahead of the phase 
transition bubble wall within the symmetric phase. In particular, we generalize ex- 
isting approaches to distinguish between chemical potentials of particles and their 
superpartners. This allows us to test the assumption of superequilibrium (equal 
chemical potentials for particles and sparticles) that has usually been made in ear- 
lier studies. We show that in the Minimal Supersymmetric Standard Model, su- 
perequilibrium is generically maintained ~ even in the absence of fast supergauge 
interactions - due to the presence of Yukawa interactions. We provide both analytic 
arguments as well as illustrative numerical examples. We also extend the latter to 
regions where analytical approximations are not available since down-type Yukawa 
couplings or supergauge interactions only incompletely equilibrate. We further com- 
ment on cases of broken superequilibrium wherein a heavy superpartner decouples 
from the electroweak plasma, causing a kinematic bottleneck in the chain of equili- 
brating reactions. Such situations may be relevant for baryogenesis within extensions 
of the MSSM. We also provide a compendium of inputs required to characterize the 
symmetric phase transport dynamics. 
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1. INTRODUCTION 

Electroweak baryogenesis (EWB) remains one of the most attractive and testable scenar- 
ios for explaining the baryon asymmetry of the universe (BAU). The BAU is characterized 
by the baryon number-to-entropy density ratio Yb =nB/s, where hb is the baryon number 
density and s is the entropy density. The value for Yb obtained from analysis of light ele- 
ment abundances in the context of Big Bang Nucleosynthesis (BBN) is consistent with the 
value extracted by the WMAP collaboration from acoustic peaks in the cosmic microwave 
background anisotropy: 



6.7 X 10-11 < Yb < 9.2 x IQ-^^ (95% C.L.) BBN ]!} 
8.36 X 10-11 < Yb < 9.32 x 10"" (95% C.L.) CMB [1, 2]. 



As observed by Sakharov [3], a particle physics explanation for this tiny number, corre- 
sponding to roughly five percent of the cosmic energy density, requires three ingredients 
(assuming a matter-antimatter symmetric universe at the end of inflation): (1) violation of 
baryon number, B; (2) violation of both C and CP symmetry; and (3) a departure from 
thermal equilibrium, i 

In EWB, these ingredients come into play if the scalar (Higgs) sector of the theory gives 
rise to a strongly first order electroweak phase transition (EWPT) at temperatures T ~ 100 
GeV. Such a phase transition proceeds via bubble nucleation, wherein regions of broken 
electroweak symmetry emerge in a background of unbroken electroweak symmetry. C and 
CP-violating interactions of fields near the bubble wall lead to the creation of left-handed 
charge that is injected into the unbroken phase, where electroweak sphalerons convert it into 
baryon number. The bubbles expand into regions of ub 7^ 0, freezing it in because sphaleron 
transitions are quenched in the bubble interiors. A strong first order EWPT is required in 
order to sufficiently quench the sphalerons inside the bubble, thereby preventing wash out 
of the captured baryon number density. 

Although the Standard Model (SM) in principle contains all the necessary ingredients for 
EWB, the effects of SM CP-violation are too suppressed to generate sufficient left-handed 
charge during the process of electroweak symmetry-breaking. Moreover, the LEP II lower 
bound on the mass of the SM Higgs boson, rrih > 114.4 GeV is too high to allow for a strong 
first order EWPT ^. Consequently, EWB can be viable only in the presence of new physics 
at the electroweak scale. In particular, augmenting the scalar sector of the SM can lead to 
a strong first order EWPT consistent with a SM-like Higgs scalar that is heavier than the 
direct search lower bound. The possibilities for doing so encompass both supersymmetric 
and non-supersymmetric scenarios, and in either case, searches for new scalars at the Large 
Hadron Collider could provide important tests (for recent work, see Refs. 0,0,0,1,1, EH 



and references therein). Similarly, the presence of new CP- violating interactions, the effects 
of which are not suppressed by light quark Yukawa couplings and small mixing angles as it is 
the case for the CKM mechanism [12], could lead to sufficient left-handed charge generation 
during a first order EWPT. Experimental searches for the permanent electric dipole moments 
(EDMs) of the electron, neutron and neutral atoms with enhanced sensitivity could uncover 



the existence of such interactions [13|, [IJ, Il5| . In addition, even though other evidence for 



^ This criterion can be evaded if CPT invariance is broken. 
^ Numerical studies indicate that in a SM universe, electroweak symmetry breaking occurs through a smooth 
cross over rather than through a phase transition 
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CP- violation may be provided by B-physics (igI . 17| (although not for the minimal field 



content of MSSM [18|), a direct measurement of the parameters fixing the relevant CP- 



violating physics will most likely require a collider beyond LHC such as the ILC |18l . |19 
In light of these prospective experimental searches for the ingredients needed for EWB, 
it is important to refine the theoretical apparatus for relating their results to the baryon 
asymmetry. 

The left-handed charge density, that biases weak sphaleron transitions and that is there- 
fore of central relevance for the computation of the baryon asymmetry, is given by the sum 
of all charge densities of left-handed quarks and leptons of all generations, 

3 

nieft = J2'yqi + ii). (2) 

i=l 

The number densities are understood to be the sum of both isospin components and of the 
three colors for the quarks. Moreover, we use the word density as a short hand expression 
for charge number density (the zero component of the vector current, which is the difference 
of particle and antiparticle number densities), and denote the densities by the symbols 
that also represent the particular particles. As outlined above, the left handed density 
gets converted into a baryon density ns through weak sphaleron transitions. The following 
formula describes baryon generation and washout ahead of the bubble wall 20| : 





ub = -3— / dz nic{t{z)e^^\ (3) 

J 

— CO 

where f ^ is the bubble wall velocity and z is the spatial coordinate perpendicular to the wall 
in the frame where the wall is at rest. Negative values of z correspond to the symmetric 
electroweak phase (bubble exterior), positive values to the broken phase (bubble interior). 
Because the weak sphaleron rate, r„s, is much slower than the rates for both the creation 
of TT-ieft and its diffusion ahead of the bubble wall, application of Eq. ([3]) is usually decoupled 
from the network of diffusion equations, a simplification that we also adapt here. Eq. ([3]) 
underlines the essential need of accurate theoretical methods of determining riieft in order to 
make quantitative predictions for Yb- 

The importance of diffusion for EWB has been emphasized in Refs. (20l . 21 . 2^ . Due to 



scatterings with the thermal bath, the CP- violating density riieft is not only localized at the 
bubble wall, but it is also transported to the region ahead of the bubble wall. Therefore, there 
remains a larger amount of time for weak sphaleron processes to turn nieft into the baryon 
asymmetry, before the it is captured by the bubble inside of which sphaleron transitions are 
quenched. 

The purpose of the present paper is twofold: First, we present a more detailed discussion 
on the derivation of the diffusion equations and the computation of the interaction rates that 
enter these. This supplements our recent publications |23l . H]. Second, we extend the net- 



work of diffusion equations to distinguish between particle and sparticle chemical potentials. 



In earlier publications (see Refs. [20|, l23|, |2J] and references therein), it has been assumed 



that the chemical potentials for particles and their superpartners are identical, a situation 
that we refer to as superequilibrium. Here, we provide numerical evidence that superequi- 
librium holds in most regions of parameter space. On the other hand, the generalization 
presented here also allows for a computation of the baryon asymmetry in parametric regions 
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where superequilibrium does not hold. As for the source of CP- violation, although we fo- 
cus in this paper on the MSSM as an illustrative case, our methods can be applied to any 
super symmetric scenario [e.g., the "next-to-minimal" Supersymmetric Standard Model). 



1.1. Existing approaches to diffusion 

In order to compute the left-handed charge, we need to derive and solve a coupled set 
of transport equations for the densities of particles that couple directly or indirectly to the 



CP- violating sources [20|, [21 



d,ji^ = -J2^rsns + Sf^ , (4) 

s 

where is the current density for particle species r, Trs are transport coefficients that 
couple the evolution of species r to the number densities of other species s, and Sf^ IS 
a CP- violating source term for the species r. 

In earlier treatments on diffusion for EWB, it is usually assumed that the only relevant 
Yukawa coupling is the one of the top quark, as it is much larger than Yukawa couplings 
of the first generations but also much larger than third-generation couplings of down-type 
quarks and leptons [20!]. While we agree with the general framework for the diffusion equa- 
tions and the strategies for analytical solutions that is described Ref. [lO] and followed in 



most subsequent work, in two recent publications [23|, |2J] , we have shown that the ratio be- 
tween the Yukawa couplings does not directly answer the question of their relevance. Rather, 
the timescale associated with the interactions induced by these Yukawa couplings has to be 
compared to the inverse diffusion length F^j^ that is characteristic for the EWPT. In super- 
symmetric scenarios, the Yukawa couplings of down-type fermions and their superpartners 
grows with tan (3, thereby enhancing the equilibration rate for these reactions relative to dif- 
fusion. We find that for tan /3 > 5 (tan (3 > 15) interactions between bottom- (r-) particles 
and the Higgs sector can in general equilibrate on diffusion time-scales. This observation 
induces important qualitative and quantitative changes to the description of the diffusion 
process in EWBG: 

• When bottom quark Yukawa couplings are in equilibrium, no net chemical potential 
associated with the axial charge of left handed third generation fermions arises. As 
a consequence, the production of densities of first generation quark densities through 
strong sphaleron processes (thermal SU(3) instantons) is suppressed. 

• The sign of the baryon asymmetry depends on the sparticle mass spectrum and on 
tan/3, and it is therefore not uniquely given in terms of the CP-violating phase. In 
particular, in parametric regions where r- Yukawa couplings are negligible but bottom- 
Yukawa interactions equilibrate, the sign changes according to whether the right- 
handed sbottom-particle is heavier than the right handed stop or not. 

• When also r- Yukawa interactions are in equilibrium, there are substantial contribu- 
tions from third generation leptons to nieft. 
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1.2. Supergauge interactions and diffusion 

The particular step towards a complete treatment of the diffusion dynamics in the sym- 
metric phase that we take in this paper is to generalize the diffusion equations to allow 
for different chemical potentials for particles and their superpartners, thereby accounting 
for possible deviations from superequilibrium. This is of importance because in supersym- 
metric scenarios the CP- violating sources generally involve interactions between supersym- 
metric particles and the space-time varying Higgs vacuum expectation values, leading to 
non-vanishing densities for the superpartners. Since the electroweak sphalerons feed on a 
net left-handed charge for quarks and leptons, supersymmetric interactions must efficiently 
transfer the non-vanishing superpartner densities into an asymmetry involving left-handed 
SM fermions. 

In the MSSM, the most important source of left-handed charge is CP-violation in the 
Higgsino-gaugino sector that gives rise to a non-vanishing Higgsino density (see, e.g. 



Refs. 25|, l26|, |27] and references therein). The source requires a non- vanishing phase 
between the supersymmetric parameter and the SUSY-breaking gaugino mass parameters, 
Mi^2- Its magnitude is largest when the difference between /i and either Mi or M2 is small 
compared to their magnitudes, leading to so-called resonant electroweak baryogenesis. In 
previous work, it is usually assumed that supergauge interactions , such as if 
are sufficiently fast, that once generates a non- vanishing Higgsino density, the latter 
immediately converts into a non-vanishing Higgs boson density. Under this assumption 
of gaugino-mediated superequilibrium, one may work with a total density H for the Higgs 
bosons and their superpartners. Yukawa interactions convert the H density into that for SM 
quarks and leptons, which are again assumed to be in equilibrium with their superpartners. 

These assumptions are indeed well justified when the supergauge interactions are in equi- 
librium. By not distinguishing between chemical potentials of particles and their super- 
partners, superequilibrium has been implicitly imposed in producing the numerical results 
presented in Refs. j23|, l^]. In this paper, we do distinguish between particle and sparticle 
chemical potentials and take accurate account of the finite interaction rates that tend to 
establish superequilibrium. In doing so, we show that: 

• Superequilibrium holds in most of the relevant MSSM parameter space. We produce 
both analytic arguments to illustrate the reasons why and numerical studies for pa- 
rameter space regions where the analytic arguments break down. 

• Superequilibrium yet may be maintained when supergauge interactions are slow, ei- 
ther because the gauginos are heavy and decouple from the plasma or because the 
corresponding three-body interactions are kinematically forbidden. This preservation 
of superequilibrium arises through a chain of reactions involving Yukawa interactions. 

• In the gaugino decoupling regime, the chain of Yukawa reactions may be broken due 
to the presence of addtional heavy (s)particles, leading to a departure from superequi- 
librium. The assumption of superequilibrium in this case may lead to an unrealistic 
prediction for Yg- Such scenarios may be relevant in extensions of the MSSM that do 
not require light gauginos for the existence of significant CP-violating sources in the 
transport equations (jl]). 
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1.3. Outline of this paper 

The plan of this paper is as follows: Section [2l culminates in the full network of Boltzmann 
equations that describe diffusion, which we present in Section [2^ Leading to this, we discuss 
how these equations may be derived within the closed time path formalism (Section 12.11) . 
list the relevant interactions within the MSSM (Section 12.21) . introduce the fully thermally 
averaged three-body supergauge interactions (Section 12.31) . discuss the thermally averaged 
Yukawa and triscalar interactions (Section l2.4l) and the particular source and relaxation rates 
for the asymmetry (Section [2.6p . Additional inputs needed are thermal masses (Section [2. 7p 
and diffusion constants (Section 12. 8p . Some simplifications in the approximation of thermal 
effects in the averaged interaction rates are discussed and justified in Section I2.10[ 

In Section [3l we present approximate analytical solutions to the Boltzmann equations. 



While in Section 13. 2[ a brief review of the discussion in Refs. [23|, l2J] is provided, in Sec- 
tion 13.11 we go beyond that and show how superequilibrium can be maintained even in 
case when supergauge interactions are quenched (e.g. through large gaugino masses or for 
kinematic reasons). 

In Section |U we provide the numerical evidence for the preceding discussions. An illus- 
trative point in parameter space is presented, where the analytical approximation is justified 
and yields reasonably accurate predictions for the densities ahead of the wall (Section l4.ip . 
By variation of tan /?, the effect of changing the strength of down-type Yukawa couplings is in- 
vestigated in Section |l]2l We refer the reader to Fig. H] which illustrate this tan /3-dependence 
on the relationships between chemical potentials for particles and their superpartners, and 
to Fig. which shows the corresponding impact on the relationships between left- and 
right-handed Standard Model fermion chemical potentials. These figures also illustrate the 
impact of Yukawa-induced superequilibrium that would persist if the supergauge interaction 
rates were set to zero. How superequilibrium is broken and how it can be maintained in the 
absence of supergauge interactions is exemplified in Section 14. 3[ Fig. [6] summarizes the final 
impact on the baryon asymmetry of several of these features. There we give Yg as a func- 
tion of tan/3 that results from the full computation with our benchmark input parameters 
and compare to the results that would have been obtained had we neglected the presence 
of supergauge interactions, the third generation lepton contributions, or departures from 
superequilibrium that arise at large tan (3. Conclusions are drawn in Section O 



2. DIFFUSION TRANSPORT EQUATIONS 

In this Section, we discuss in detail the derivation of diffusion transport equations. In 
comparison to earlier treatments, we generahze these equations to distinguish between par- 
ticle densities and the densities of their superpartners. 



2.1. Three-body rates: general formalism 

We derive the diffusion transport equations for EWB using the closed time path (CTP) 
Schwinger-Dyson equations. Although one may use conventional kinetic theory for this 
purpose, we adopt the CTP framework as it allows one to systematically include higher- 
order corrections and the effects associated with departure from adiabatic quantum evolution 
(for a detailed review of the CTP framework as applied to EWB, see our earlier work in 
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Refs. 25|, l28|). The CTP transport equations for bosons and fermions have the form 26l. |29|: 

dnB 



dXo 



(X) + V-jB(X) 







^(X) + V-j^(X) = -Jd'z r dzn Tr 



S> (X, z)G<{z, X) - G>iX, z)S< X) 
+G<{X, z)J:>{z, X) - S<(X, z)G>{z, X) 

E>(X, z)S<iz, X) - ^>(X, z)S<(z, X) 
+^<(X, z)E>{z, X) - S<(X, z)S>{z, X) 



(5) 



(6) 



where = {jiB^jB) and = (ni?, jV) are the boson and fermion current densities, respec- 
tively. The functions G^' ^{x,y) and S*^' ^{x,y) are elements of the 2x2 matrix of CTP 
boson and fermion propagators: 



S^\x,y) = {Tr (y) ) 



(7) 
(8) 



where (...) denotes an average over the physical state of the system, T-p is a path ordering 
operator, and the indices a, b denote the branch of a closed time integration path running 
from — oo to +oo (the "+" branch) and back to — oo (the "— " branch) the fields inhabit. In 
the case of the bosonic Greens functions, one has 



= G\x,y) = {T [<P{x)<p\y)]) 

G+-{x,y) ^ G<(x,i/) = (0t(^)0(x)), 

= G>(x,i/) = (0(x)0t(^)), 

G-{x,y) ^ G\x,y) = {f [<l>{x)<p\y)]) 



(9) 
(10) 

(11) 
(12) 



while the corresponding expressions for the elements of S°'^{x,y) contain the appropriate 
factors of —1 to account for fermion ant i- commutation relations. The self energy functions 
S^' ^{x, y) and S^' ^(a;, y) give the corresponding one particle irreducible corrections to the 
free inverse CTP propagators. 

The source terms on the RHS of Eqns. ([5]) and ([6]) can be obtained by computing 



x, y) and S^,' (x, y) order-by-order in perturbation theory. In general, doing so re- 



quires knowledge of the non-equilibrium distribution functions. However, the presence of 
a hierarchy of scales allows one to expand these functions about their equilibrium values 
in powers of appropriate scale ratios. As discussed in Refs. [25|, l28j, these scales include 



a decoherence time, r^, associated with the departure from adiabatic dynamics; a plasma 
time, Tp, associated with mixing between degenerate states in the plasma; and an intrinsic 
quasiparticle evolution time. Tint, associated with the time evolution of a state of a given 
energy. For the dynamics of the electroweak plasma, one finds that Tint <^ Tp <^ r^, leading 
to a natural expansion in the scale ratios: Sd = Tint/r^ ~ f^fces/i^ and Ep = Tint/Tp ~ Tp/uj 
with Vw being the bubble wall expansion velocity, k~g being an effective length scale (such 
as the wall thickness), Tp being a thermal quasiparticle "damping rate" in the plasma, and 
u being the quasiparticle frequency. Here, we include in the damping rate any an process 
involving emission and absorption from the thermal bath. In addition, the plasma is rela- 
tively dilute, so that the ratio of chemical potentials to temperature, = fi/T, provides an 
additional expansion parameter. 
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In terms of the e parameters, the leading contributions to the RHS of Eqns. ([5]) and ([6]) 
occur at O^e^). Specifically, the CP- violating sources include effects of order x Sp or 
e% while the CP-conserving sources arise at order SpX e^. In this context, the supergauge 
interactions generate terms of the latter type. As discussed in Refs. jH, H], the presence 



of the scale hierarchies embodied in the e parameters allows us to adopt the quasiparticle 
ansatz for the CTP Green functions and to work near chemical and kinetic equilibrium when 
computing these terms. To this end, we compute the the terms on the right hand side of 



the transport equations fl5|6l) following the procedure used in Ref. [28| for the calculation of 
the Fy-type terms. The bosonic CTP Green functions entering the computation are given 
by 

/d^k ( \ 
j—^^ \\ + /ij(/i;o, /ii) j a(A;o, k) , (13) 

/i?(^o, /ii) Pi{^Q, k) , (14) 

with spectral functions 
pIU, k) = i/^k - ^k) - i[W + i^k)] (15) 

= ^[f i i \-{ 1 i 

2cjfc \ki^- ujk^ie k^ + ujk + ie J yk^ - ujk - ie ko + Uk - ie 



where = y |k|^ + m^. The pj(fco,k) can be appropriately modified to take into account 
collision-broadening and thermal masses. Using the expansion in e parameters introduced 
above, it suffices to take the distribution functions to be close to the equilibrium form 

fB{ko,Pi) = nsiko^Hi) + 0{ed/ep) , (16) 

where nsiko, fJ^i) = l/[e*^'^"~^'^/^ — 1] and fii is a local chemical potential . Here we neglect 
the terms of order e^/Sp ~ v^k^s/Tp <^ 1. 

Similar expressions are obtained for the fermion Green functions: 

S>{x,y) = mx)i,{y)), (17) 

S<{x,y) = -{^{y)^{x)), (18) 

which can be expressed as 

S\x, y) = j ^^e-^H^-y^ f^{ko, ^)p{ko, k) (^ + m) , (19) 

where A denotes either ">" or "<" and and with the functions 

f^{ko,n) = l-nFiko-n), (20a) 
fP{ko,fi) = -nF{ko-fi), (20b) 

with nF{ko,fii) = l/[e^''o-i...)/T _^ 
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2.2. MSSM interaction Lagrangian 

In order to calculate the transport coefficients in the Boltzmann equations, we ffist iden- 
tify the relevant interactions in the MSSM Lagrangian. These interactions, denoted by /^mt, 
can be divided into three classes: 

Bilinear interactions that arise when the neutral Higgs bosons acquire vacuum expectation 
values (vevs) are 

'Cm = -ytt*ji ti {At f „ + ^i*Vd) - yt Vu in Pl ti (22) 

- Vb h*R bL {Ab Vd + n) - Vb Vd Ir Pl &l 

- Vt tl {A^ Vd + fJ.*Vu) - Vr Vd tr Pl tl 

- ^ ^^o{vd Pl - e'^"\^ Pr)^^ - ^ ^ j^oivd Pl + e"^"\u Pr)'^wo 



- g2 -^H+^Vd Pl + e'^^ + h.c. 



These terms, which result in squarks, quarks, sleptons, leptons and Higgsinos scattering due 
to the spacetime-dependent Higgs vevs Vu{x) and Vd{x)-, contribute to CP- violating sources 
and CP -conserving relaxation rates Ym and Yu- The calculation of has received 
the most attention, both in the CTP approach and in other frameworks; however, there 
remains a significant dispersion in the recent calculations of separate groups I2I, i^, Q H 



331]. The CP-conserving relaxation rates have been estimated in Ref. [20[, and rigorously 
calculated and studied with CTP methods in Ref. [25| in a manner consistent with the 
computation of the CP- violating sources. 

Trilinear interactions proportional to the top Yukawa coupling yt are 



(23) 



ICy. =-ytt\ tL {At Hi + /x*i70*) + yt bL {At H+ - ^i*H,*) 
+ yt {H+ tR Pl bL - Ir Pl II) 
+ yt e'^- {tR h Pr + tR II Pr ^1+) 

+ yt e-'^^ {tL tR Pl ^ ho +bL tR Pl + h.c. . 

For the bottom Yukawa coupling y^ the corresponding interactions are 

^y. = - ybb*RbL (A i/o + ^^*Hl*) + yj^f^ (A - (24) 
+ yt {H^ Ir Pl tL - H^d bR Pl bL) 
+ Vb (-^R h Pr^ho+ bR h Pr^h+) 
+ yb (bi bR Pl - tl bR Pl ^ |+) + h.c. , 



where we have employed notations and rephasing conventions according to Ref. [25|. The 
corresponding interactions Cy^ for third generation (s) leptons follow when replacing bR — >■ tr 
and tL ^ i^T 'Cyi^. 
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In general, all interactions mediated by third generation Yukawa and triscalar couplings 
may have a sizeable impact on the result of EWB, such that we consider 



Vb 



■'Vr 



(25) 



These terms, which are Yukawa interactions and their supersymmetric counterparts, lead 
to transport coefficients generically denoted Fy. These equilibration rates are potentially 
important in that they communicate CP-asymmetries from the Higgs sector to the quark sec- 
tor, biasing sphalerons to produce a baryon asymmetry. The dominant absorption/emission 
contribution to Fy has been studied within the CTP framework in Ref. 28| : the sub- 



dominant scattering contribution to Fy has been partially calculated in Refs. [20|, l21 
Supergauge interactions are 



C 



V 



_ 9i_ 

_ 92_ 

-92 
_ 92_ 

V2 
_ _9i_ 

3^2 

-g3V2 



^n^-H,* Pl + e^*"'if+ Pr)^^o + ^ m^HT Pl + e'^"^' Hi P«)vI/^„ 
^n^{H'd Pl + e''*''^"Hl Pr)^^^ + ^y^^{H-* Pl - e^^"'i/+ P^)vl>|„ 
u\* Pl ul - d'* Pl d\, + uf, ^^,oPl z^l - ^r. ^woPl e ' 



(26) 



'wo ^L «L 



H?o^L l^L - ^WO-^L fiL 



^l^^PLV^L + d'l^nPLd 



B 



+ ^ [K^nPLi^l + ^l^-BPLel] 



u'l X'^m, Pl u\ + dl A'^ ^r?, d\ 



g, V*^ 4* Pl u\ - V,, Pl di - g, ^l^^, Pl i^l - g^ vt^^ Pl e\ 



u], r^^^P^vl/^ + d^n >^''d^RPL^l 

d^nd'nPL^B 



+ ^g, [u),n^^P^^^]-y^g, 



+ V2gi [?^4Pi^5] + h.c. 



where A° and A" are the generators of 3 and 3 of SU(3). These supergauge interactions - 
the supersymmetric version of gauge interactions in the SM - lead to transport coefficients, 
generically denoted by F^, which tend to equilibrate the chemical potentials for particles 
and their superpartners. All previous studies have assumed the limit Fy — > oo, which leads 
to superequilibrium. One of the purposes of the present work is to calculate Fy and to solve 
the Boltzmann equations without the assumption of superequilibrium. 

In employing these interactions to compute the supergauge equilibration rates, we work 
with the mass eigenstates of the unbroken phase: gauginos and Higgsinos (rather than 
charginos and neutralinos) , left- and right-handed quarks and squarks, and Higgs scalars. 
Deep inside the bubble, this choice is clearly not appropriate, owing to large flavor mixing 
induced by the non-zero Higgs vevs. A proper treatment of this flavor mixing requires 
a modification of the transport equations that allows for an all-orders summation of the 



spacetime varying Higgs vevs [33|, |3J|. In the absence of such a treatment, we will work in 
the "vev insertion" approximation, defined by assuming that the dynamics of chiral charge 
production are dominated by the region near the phase boundary and that the Higgs vevs 
in this region are small compared to the temperature and slowly varying (i.e., the wall is 
relatively thick). Under these assumptions one may treat the flavor mixing perturbatively. 
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We find below that the particle densities are generally largest in magnitude near the phase 
boundary as one would expect if the vev insertion approximation is valid. Nonetheless, we 
emphasize that our specific numerical conclusions are provisional and await a more complete 
treatment of the flavor mixing dynamics in the broken phase within the bubble. 

Within the spectrum of unbroken phase eigenstates, the Higgs scalars provide additional 
complications. To appreciate the difficulties more clearly, suppose for simplicity the light- 
est Higgs mass eigenvector field is a fixed linear combination of Hu and Hd as we pass 
through the bubble wall. Then the naive dispersion relationship for modes in the WKB 
approximation will be 

= ^/W+V^ (27) 

where V{(j)) is the thermal effective potential and v{x) = {(f)) solves the Euler-Lagrange equa- 
tions involving V'{(f)) leading to the cancellation of tadpoles. For states with sufficiently large 
momentum, the fact that V"{(j)) < in a particular (f) range will be unimportant. However, 
sufficiently soft modes will become unstable if V"{(f)) < and may lead to significant back- 
reaction corrections to v{x) on time scales of order Imco^^ and can also lead to particle 
production effects. Note that this statement of the problem is a bit more subtle than it 
seems because v{x) is not spatially homogeneous while the thermal effective potential V{v) 
obtained by the usual construction methods corresponds to the energy of the system with v 
fixed to be spatially homogeneous (e.g. see (35|). 

In this paper, we will simply settle with an estimate of the error incurred by neglecting 
this inhomogeneous bubble profile effect on the Higgs transport equations, and leave a more 
detailed analysis to a future work.^ If we denote the largest magnitude of V"{(j)) < region 
when the temperature is near the critical temperature as kc = max[|V""(0)|^/^] , the effect 
of neglecting of the V"{(j)) < region on the Higgs transport equations can be estimated by 
the following fractional thermal distribution number density: 

n H{oo,m/Tc) 

where 

H{kjT,m/T) = / dxx^/x^ - m'^/T'^ [exp {x) - 1]'^ (29) 

Jm/T 

and where m is the unbroken phase mass. The function Uc/n is a monotonically increasing 
function with kc/T and a monotonically decreasing function with increasing m. Although 
the exact value of kc/T^ is model dependent, if we take 1 ^ kc/Tc ^ 2 and take m/Tc = 0.7, 
we arrive at fractional correction estimate of 

Ti 

0.09 < — < 0.35. (30) 
n 

Hence, we can expect order 10% of the Higgs scalar density will have significantly distorted 
phase space distribution due to the instability of these modes. Although not completely 
negligible numerically, these effects can be seen as refinements on the order unity effects 
presented in this paper. 



^ Any large effects coming from this is unlikely to be computable analytically. 



(28) 
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In the remainder of the paper, we will make the simplifying assumption that quadratic 
fluctuations about the classical solution to the field equations can be approximated by a set 
of scalar fields (-ff^, , H^) with mass terms 



H7' 



(31) 



and the same for H^^) but with b + 6b 



— (6 + 6b)- The terms 6u,d,b denote finite 
temperature corrections which lead to a minimum in the Higgs scalar potential at f « = f d = 
0. We can re-express this mass matrix using the minimization conditions for electroweak 



symmetry breaking at T = [36 



m 



Hu 



+ l/ip = cos^ /5o + -m| cos2/?( 



m 



2 

sin^ — 2^^! ^'^^ 2/5o . 



h = rrij^ sin cos , 
where and rriA are the Z and pseudoscalar Higgs boson masses at T = and 



tan/3o = — 

Vd 



(32) 

(33) 
(34) 

(35) 



T=0 



Therefore, the mass term becomes 



^ / +1 N / cos^ /5o + cos 2/?o + ^« j 
[^u,^d)y lmlsm2po + 6b m^sin^ 



im^ sin 2/?o + ^6 



sin Po — |m| cos 2/3o + 6d 



H-d' 
(36) 

The eigenvalues, corresponding to the charged Higgs scalar masses in the unbroken phase, 
are 



m 



Hi,. 



ml + {6u + 6d)Ty {{m\ + m|) cos2/?o + {6^ - 6^)) + (m^ sin2/?o + 25^)^ 



(37) 

We see that if 6u,d,b 0, then m|^^ < — a consequence of the fact that electroweak 
symmetry is broken at low temperatures. Furthermore, the mixing angle a, defined by 



Hi 



COS a sm a 
- sin a cos a 



is given by 



tan 2a 



m\ sin 2/5o + 26b 



(38) 



(39) 



cos 2/3o + ^1 cos 2/5o + 6u — 6d 

(The neutral Higgs mass matrix differs from the charged Higgs mass matrix only hy h + 6b — *■ 
— (6 + 6b)'., the mass eigenvalues are the same and the mixing angle differs by an overall 
sign). We emphasize that we are diagonalizing the Higgs potential about its minimum 
at {Hu) = (Hd) = in the unbroken phase, as opposed to the usual zero-temperature 
treatment [36| ; our results for rn^j^ ^ and a do not simplify to the zero-temperature Higgs 
masses and mixing angles in the limit 6u,d,b 0. 
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If we assume that m^, 5u,d,b, and tan/^o ^ 1; then 



' - Su- \ml (40) 
"^^2 - m\ + + 5rf (41) 

and a ~ 1/ tan/3o. In our analysis, we will use Eq. (!39ll . In general, this mixing angle will 
be spacetime-dependent due to the appearance of terms in the mass matrix proportional to 
Vu{x), Vd{x). However, it has been found that 

A/3 ^ /3iT) - /3{T) (42) 

z — ^oo z — i-— oo 

is numerically small: A(3 ^ 10^^ [37]. Therefore, we neglect the spacetime-dependence of 
the Higgs mixing angle. (In extended supersymmetric EWB scenarios, it is conceivable that 
A(3 might be larger, necessitating a proper treatment of spacetime-dependent Higgs mixing.) 

The purpose of the preceding analysis is to motivate realistic masses and mixing angles 
in the Higgs sector. We defer a rigorous determination of Higgs boson masses and mixing 
during EWB to a future study. 



2.3. Supergauge equilibration rates 



Supergauge interactions generate three-body absorption/decay processes as illustrated in 
FIG. lU^a). These processes drive the plasma toward superequilibrium, the condition where 
the chemical potentials for a particle and its superpartner are equal. Following closely the 



derivation of Fy that is presented in Ref. [28|, we calculate the fully thermally averaged 



rates Ty in the on-shell limit. 




FIG. 1: Examples of absorption/decay (a) and scattering (b) processes which lead to superequi- 
librium. 

We compute the supergauge interaction rates F^ arising from emission/absorption pro- 
cesses in the thermal plasma [FIG. [U^a)]. Each supergauge interaction term in Eq. (l26l) can 
be cast in the general form 

Cint = <p^{gLPL + 9R Pr) V + h.c. (43) 

for gaugino V, and (bosonic, fermionic) superpartners {4>,ip)- These interactions give con- 
tributions to the RHS of Eqns. (!5|6|) of the form 

d,f4X) = -d,j;{X) = Sy, (44) 
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where 



(45) 



The functions Ip and Tp are defined to be 2 



Xpimi, m^, ma) = -^^^3^ ("^1 + ml - m^) / dui j du^ (46) 
X |nB(u;<^) [1 - nF(u;i)] np^ui - u;,/,) [6'(mi - m2 - m^) - 6{m^ - mi - ma)] 
+ nB{uJcj,)nF{^^i) [1 - ^F(t^i + cj0)]6'(m2 - mi - m<^) 
with integration hmits on given by 



and lastly 



2m^ I ^' - ^ ^' (47) 
± [{ul - ml) [ml - (m2 + m^)^) (m^ - (ma - m^)^)] ^''^j; 



lF{mi,m^, ma) = 5 2 2-^(^1, m^, ma) . (48) 

mi + mg — mi 



The gaugino chemical potential fiy only appears on the RHS of Eq. (H^ when the gaugino 

is a Dirac fermion (i.e., for W^). For Majorana gauginos, no gaugino chemical potential ap- 
pears. Although a Majorana chemical potential does not correspond to a conserved Noether 
current, it can arise as a deviation from a pure Fermi-Dirac distribution when annihilation 
processes are out of equilibrium. However, such a deviation does not contribute to Eq. (l45l) 
due to CP-symmetry. The charge current densities and corresponding chemical potentials 
for Dirac fermions and complex scalars are all odd under CP. However, a Majorana chemical 
potential is even under CP. Therefore, in the limit that we can neglect CP-violating phases 
in our interaction rates, Majorana chemical potentials do not contribute to the Boltzmann 
equations for charge current densities. We have explicitly verified that a non-vanishing Ma- 
jorana chemical potential ultimately cancels from Sy to linear order in and at zeroth 
order in CP-violating phases (pcp- Physically speaking, an excess of Majorana gauginos V 
does not bias a charge asymmetry in ip and since the rates for V <-> and V ^ ipcf)^ are 
equal"^. 

At present, we work exclusively within the MSSM, where existing measurements constrain 
(pcp ^ 1; consequently, we neglect these contributions. In extensions of the MSSM where 
additional CP-violating phases are less constrained and may be large, the non-equilibrium 
dynamics of Majorana fermions may be important. We also emphasize that these arguments 



^ In leptogenesis scenarios, a heavy Majorana neutrino can bias a chiral lepton asymmetry, as long as the 
relevant CP-violating phases are non-zero. This is consistent with the statement that, in the present 
discussion, transport coefRcients which couple Majorana chemical potential /iy to /i^, /z^ arise at order 
©(sin (j)cp)- 
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apply to the degrees of freedom that are present in the symmetric phase. This is consistent 
within the framework of the present paper, as our main focus is the diffusion process in the 
symmetric phase and as we calculate the source and relaxation terms in a pertubative mass- 
insertion scheme. Note that in the broken phase, the Higgsinos, which are treated as charged 
particles in the symmetric phase, mix e.g. with the Binos, which are Majorana particles. 
The resulting neutralinos are Majorana particles, but due to their Higgsino component, their 
out-of-equilibrium dynamics is crucial for EWB. A non-pertubative treatment of the mixing 
in the broken phase will be subject of future investigations. 

We can relate the chemical potentials to charge number densities n via 



^ - / i-T-^PT^ - ...,J.,r , . 1 = ? k{mlT) /. + 0{el) , (49) 



(27r)3 \e('^p"^)/'^ ± 1 e('^f +^)/^ ± 1/ 6 



where k{m/T) is a statistical factor [25 



k{m/T) = A;(0) / dxx -^r J - w? /T^ (50) 

^ Jm/T [e'' ± 1)^ 

with k{0) = 1 for chiral fermions, A;(0) = 2 for Dirac fermions and complex scalars, Cf{b) = 
6(3) for fermions (bosons), and the + (— ) sign for fermions (bosons). Therefore, we can 
write 

for Majorana V, where 



V J"2 



6 I I I I 

( I + Is-rI ) ("^V''"^</''"^y) +2 Re (fifLS-^j) (m^,m<^,m^) . (52) 



If V is Dirac, as is the case for W^, the situation is more subtle. Instead of Eq. (IFIil . 
when inserting Eq. ( |52l) into Eq. ( H5l) . we have 



S = rl^'*) (^_^_^\ , (53) 



^ ^ kfj) ky 

meaning that that non-equilibrium dynamics of does not in general decouple from the 
dynamics that produces Yb. Let us consider the contributions from Syj^± to the Boltzmann 
equations (jS]) for the third generation LH quarks: 



Now we define 





















l{^'h-^^b,) 



(55a) 
(55b) 
(55c) 
(55d) 



17 



With these definitions, we obtain from Eqns. ([5 



X (A/iq + A/ig-2/i^±) 



(56) 



and 



Ncgl \Vt 



tb\ 



If {mt^ ,mi^,m^±) - If {mb^ , mj^ , ) (/i, - /i,) 



L-l'^tL-l^W±> 



(57) 



X 



In the unbroken phase, where Vu = = 0, we have mt^{T) = rrib^iT) = niq and m^^{T) 
(T) = mg-, so that the differences 



If {mt^ ,m~^^,m^±) - If {rrib^ , m^^ , 

while the sums 

{nit^ ,m^^,m^±) +If {nib^ , m^-^ , 



2lF{mg,mq,Tn^±) . 



(58) 

(59) 



Therefore, to the extent that we can neglect isospin-violating mass differences within the 
bubble, isovector asymmetries such as A/Xg decouple from isoscalar densities such as fig. Al- 
though we have chosen only only one interaction as an illustration, we have verified that the 
decoupling of isoscalar and isovector densities occurs for all Yukawa and supergauge inter- 
actions (l23l[26l) . Since Yb is produced by weak sphalerons sourced by the chiral asymmetry 



^left 



i=l 



(60) 



itself an iso-scalar density, isospin-violating asymmetries decouple from the determination 
Yb- In particular, only couples to isospin-violating asymmetries and will decouple from 
the dynamics of nieft. Therefore, we consider the closed set of Boltzmann equations for the 
following twenty-five densities: 





H = 




(61a) 


qi,2 = n„i,2 + n^i,2 , 






(61b) 


q = q3 = nt^+ , 


Q = 




(61c) 


«i,2 = n^i^ , 






(61d) 


t = Us = rif^ , 


t = 




(61e) 


0^1,2 = n^ia , 






(61f) 


b = d3 = nbj^, 


b = 


d2 = n~^^ . 


(61g) 


£ = i3 = nr^+ n^r , 


1= 




(61h) 


r = n^j^ , 


T = 


^rn ■ 


(61i) 
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We defer to future work the study of isospin breaking effects within the bubble wall which 
couple isospin-violating asymmetries to those listed above. Note that within the context 
of non-supersymmetric models, a study of deviation from isospin equilibrium has been per- 
formed and the effect has been found to be quantitatively relevant 38[ . 

The supergauge equilibration rates that enter the Boltzmann equations for these densities 
fl6Tl) are 



p(52,H) 



-p{u,u) 
V 

V 

V 



2^2 



lFimg,mH^,mg) - sin 2a Jir(mjj, mj^^, m^) 



+ 



2^2 



lF{mg,mHi,my^^) + sin 2a Jir(m^, m^i, m^' 



2^2 



lFim^,mH2,'mj^) + sin 2a JF(m^, m^a, "^^^ 



+ 



2^2 

_ 2 Nc gl , 
~ 32^2 -LF[mg,mg,m^ 



18 Ncgl 



2"2 



mg,mg, 



) 



' 2^2 

ANcgf 



-1f{ 



•qi '"'qi '"'gJ 1 



a) 



3T2 



If {mu.rriu.m^) + 

If {md,m^,m^) + 
ISgl 



2^2 

Q{N'c - 1) 9l 



If [mu, rriu, mg) 



2"2 



If {rrid, m^, uiq) , 



2^2 



2^2 



(62a) 

(62b) 

(62c) 

(62d) 
(62e) 
(62f) 
(62g) 



where in Eqns. fl62b-e) we have omitted a generational index since these expressions are 
identical for all generations. 



2.4. Yukawa, tri-scalar, and Higgs vev induced interactions 

Presently, we summarize contributions to the Boltzmann equations that arise from in- 
teractions in Eqns. fl22| [23]) . Yukawa and SUSY-breaking tri-scalar interactions fl23l) lead 
to transport coefficients in the Boltzmann equations which couple Higgs, RH quark, and 
LH quark supermultiplets. We assume that the tri-scalar A-terms are proportional to the 
corresponding Yukawa coupling 36|. For example, the term 

Cy D yt e'^- {tn h Pr + tR h Pr ^ |+) (63) 

in Eq. (l23l) leads to a contribution to the Boltzmann equations for densities t, g, H of the 
form ^ 
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We now list the complete set of equilibration rates [28| arising from Cy: 

j^{t,q,Hi) '^^NcVt I • * I /I |2 -r ^ ^ 

L y = — — Isma/i + cos a At I 1b [mj,mj,mHi) , 

Ty = — ^ — \cosafi -smaAtl Is [m^,mg,mH2> , 

j.{t,q,Hi) '^2Ncyt 2 T ( \ 

YY = — Y2 — alF{mt,mg,mHj , 

iy — — — sm alp [rnt^mq^mHi) , 



^ Y 


2^2 


r(&,9,J^i) 




i y — 








L Y ~ 






-1 O A 7 2 

^2Ncyt 


L Y — 




rib,q,Hi) 


1 O A 7 2 

12 A^c yt 


L Y ~ 






l2Ncyl 


i y — 


2^2 




l2Ncyl 


i y — 


2^2 




12iVci/2 


i y — 






12iVcy2 


i y — 


2^2 




12 Ncyl 


i y — 


2^2 






i y — 


2^2 


J^{r,l,H2) 


12 Ncyl 


i y — 


2"2 




12iVci/2 


i y — 


2"2 



(65a) 
(65b) 
(65c) 
(65d) 
(65e) 
(65f) 

(65^ 

(65h) 
(65i) 
(65j) 

(65k) 
(651) 



I cos a /i* — sin a 1^ Xb {rn^f^ mj, rriHi) ) (65m) 

I sin a n* + cos aA^-flB {m^, rnj, rriH^ , (65n) 

Xp (^mg,mr,me) , (65o) 

siii^ a Xp {rrir, me, m Hi) , (65p) 

cos^ a Xp {mr,mi,mH2) , (65q) 

Xp [mr,mg,m^) , (65r) 
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where 



l^nBiujR) [l + nB{ujL)]nB{ujL - ujr) [9{mR - rriL - tjih) - 0{mL - niR - mj/)] (66) 
- nB{ujR)nB{ujL) [l + nB{uJL + ujR)\0{mH -niR- itil)^, 



with integration hmits given by 
I 

2m 



''R 



(67) 

± [{uj\ - ml) {ml - {mi + mnf) {ml - {mi - mn)^)] ^ ^ 



These interactions communicate the effects of CP-violation from the Higgs(ino) sector to the 
(s)quark sector. If supergauge interactions are assumed to be in equilibrium, only the sum 
of Fy rates listed above enters the Boltzmann equations that determine Yr- This procedure 



has been apphed in our recent pubhcations [23|, |2^. Without this assumption, we must 
distinguish between each rate. For example, the rate for the transfer of CP-violating effects 
from H t fl65cl) will be different from the rate of transfer from H ^ t (I65f|) ; this difference 
may impact Yr if transfer between t ^ t is inefficient, cf. Section H73l 



2.5. Four-body and ofF-shell contributions 

In addition, four-body scattering interactions, such as those illustrated in FIG.[T](b), will 
also contribute to the equilibration process [2l|. Although the associated rates are phase 



space suppressed and are higher order in the gauge couplings than the three-body rates, they 
can become leading order when the three-body processes are kinematically forbidden. At 



the same order in gauge coupling constants, also off-shell processes contribute [39|, |40|, |41 
which are loop instead of phase space suppressed. 

A calculation of these contributions to the interaction rates may hence be important 
for EWB in regions of parameter space where three body interactions are kinematically 
forbidden. We anticipate a full calculation that takes into account the masses in the Higgs 
sector and of the sfermions and the gauginos, which typically can be of the same order 
as the temperature, to be rather involved and do not perform it within the present work. 
For our parametric examples, we therefore avoid here parametric regions where three-body 
rates are kinematically forbidden as much as it is possible. The only exception that we 
make concerns the interaction between Hi, t and q, where it is hard to avoid kinematic 
suppression for three-body interactions and to allow for electroweak symmetry breaking at 
the same time. We specify the estimate for the four-body interaction that we adopt for that 
case in Section 14. 1[ 



2.6. Sources and relaxation terms 



CP- violating sources can arise from the Lagrangian terms fl22l) involving squarks and 
Higgsinos. Physically, they correspond to CP-asymmetric reflection and transmission rates 



21 



for squarks and Higgsinos scattering from the bubble wall. In Ref . 2^ , only contributions to 
source terms that become resonant are calculated; off-resonance these contributions become 
comparable to other terms that have not been calculated. Therefore, we do not include 
in the Boltzmann equations the chiral relaxation rates for squarks and Higgs bosons which 
are non-resonant for the present choices of parameters, deferring their inclusion until a 
complete calculation of these rates has been accomplished. We comment on the validity of 
this approximation in Sec. HI In the present work, we assume a resonant Higgsino source 
with = Ml = 200 GeV, which is sufficient to produce the observed BAU. The precise 
formula for is given in Ref. 25 . 

The choice of Higgsino and Bino as the source of resonant CP- violation rather than of 
Higgsino and Wino is motivated by the recent evaluation of the two loop electric dipole 
moments (EDMs) of the electron and the neutron in the limit of heavy first two generation 



squarks and sleptons [42|. The results indicate that the size of the CP- violating phase 



between /i and Mi is much less constrained by the non-observation of permanent EDMs 



than the phase between ^ and M2. Therefore, Bino-driven EWB [23|, |2J, |43| may prove to 
be a more viable option, particularly as the sensitivity of EDM searches improve. 

In addition to the source terms, there exist chiral relaxation rates that also arise from 
particles scattering with the Higgs background field fl22l) and that tend to wash out CP- 
violating asymmetries. These processes exist in both the quark and Higgs sectors and are 
generically denoted as Vm and Vh, respectively. In the present work, since we must dis- 
tinguish between the densities for superpartners, we must likewise distinguish between the 
quark and squark contributions to F^, and between the Higgs boson and Higgsino contri- 
butions to Th. In the Boltzmann equations to follow, we write down terms for all relaxation 
rates that are possible within the broken phase, both resonant and non-resonant. For the nu- 
merical examples in Section H] however, we only include the contributions that may become 
resonantly enhanced: the Higgsino-Bino (Higgsino- Wino) contribution to Th, resonantly 
enhanced for |yu| = Mi = 200 GeV, (|/i| = M2); and the matter fermion contributions to Fm, 
which are resonantly enhanced since the thermal masses satisfy the relations \mq — mt\ <^ T, 
\mq — mh\ <C T and — mr\ <C T. 



To facilitate comparison with Ref. [25|, we note that we apply the slight changes of 
notation f|'^ = F;,, FM(t,g) = (6/r2)F,-, Tuit^q) = {<olT^)T~. 



2.7. Thermal masses 

The masses that appear in the expressions for the equilibration rates Fy, Fy, in the 
relaxation terms F u and F m and in the source terms are understood to include thermal 



corrections, which we take from Ref. [2J], assuming light right-handed stop, sbottom and 
stau particles. While we treat the thermal masses as Dirac mass terms for fermions, strictly 
speaking, the thermal masses characterize a modification of the dispersion relation only 
and imply no chiral symmetry breaking in the symmetric phase. We provide a justification 
of this procedure through a discussion of the numerical inaccuracies that are incurred in 
Section EH 
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2.8. Diffusion constants 

When generalizing the distribution functions fs and fp to be dependent also on the 
spatial momentum p in an anisotropic way, it is possible to derive the Fick diffusion law 

X = -DxVX (68) 

from kinetic theory, where X^^ = {X, X) denotes the number density current for any of the 
particle species considered here. The quantity Dx is the diffusion constant for the species 
X and can be calculated from the moments of scattering matrix elements of X. For our 
purposes, it is useful to recast the diffusion law (168!) by boosting (non-relativistically) to the 
frame where the bubble wall is at rest as 

d^X^ = v^d,X-Dxd^,X. (69) 



The particular values that we use here for the diffusion constants are taken from Ref. [21[ 
and are summarized in Table [T] in Section HJ where we present results from our numerical 
studies. As a simplification, common diffusion rates for particles and sparticles are taken, 
such that e.g. Dq = Dg = Dg. The physical interpretation of Eqs. ( |68l[69ll is simple: Within 
the plasma, a particle gets rescattered and undergoes a random walk. These scatterings are 
more frequent for colored particles such as quarks than for particles that interact through 
electroweak and non-gauge couplings only, such as the r-lepton. Therefore, Dj^ is larger than 
Dq, which implies that r-leptons diffuse farther away from the bubble wall than quarks. 



In Ref. [21|, the diffusion constants have been estimated using an incomplete set of 



scattering matrix elements. In particular, it has been pointed out in Ref. [39|, |40|, that 
infrared sensitive t-channel processes have been neglected. In general, one should therefore 
expect an order one inaccuracy within the results of Ref. [2l[. The results in Ref. 



apply however to gauge theories with massless fundamental fermions only, which is why they 
cannot be applied dire ctly to sfermions and to Higgs bosons. In order to quote a value for 



comparison, from Ref. [39|, |40| we infer Dg = 4.7/T for six quark flavors, while from Ref. |21| . 
we adapt Dq = Dg = Dg = Q/T. In order to achieve predictions of accuracy better than 
order unity, it will be necessary to reevaluate the diffusion constants in the future. 



2.9. Boltzmann equations 



We now present the full set of Boltzmann equations: 

^ \kt kg kui) ^ \kt kg k^j ^ ' 

\kj kg kg J ^ \kj kt J \kj kg 
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The rate for strong sphaleron transitions is given by [20|, |44, |45 
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In order to facilitate a comparison with results appearing previously in the literature, 
we give the our values for the statistical factors in the massless limit: kg = kg/2 = 6, 
kuR = kdR = kuji/2 = k^j2 = 3, kn, = kn^ = A and kg = 4. 

Before solving the system of Boltzmann equations in Sees. El IH we conclude the present 
section by considering two issues related to the equilibration rates Tyy derived above. First, 
we consider the impact on these rates from a more rigorous treatment of massless fermion 
propagators in a thermal plasma. Second, we present a numerical comparison of Ty and Ty 
rates to demonstrate that one in general does not have Ty 3> Fy. 



2.10. Thermal fermion propagators and particle/hole modes 



For three-body processes involving scalars and fermions with masses of order T, the forms 
for the Greens functions given in Eqns. fll3|14p (and the fermion analogs) suffice for the 
computation of the three-body rates. In the case of fermions that are massless (or nearly 
massless) at zero temperature, the structure of the thermal Greens functions becomes far 
more complicated. The renormalized fermionic spectral functions contain additional poles 
- so called "hole modes" - generated by mixing between single and multiparticle states in 



the thermal bath [46|, |47[ . The resulting propagators are given by [25|, 148 
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where k is the unit vector in the k direction, and 
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where iD±{ko, k, p) are contributions to the inverse, retarded propagator proportional to 
(70 T 7"k)/2 arising from interactions. The function F{ko, k,fi) gives the non-pole part of 
the propagator, and k = |k|. The mass of these particle/hole excitations is neither zero, nor 
equal to the standard thermal mass (denoted m-r) ; rather it is a /c-dependent function given 
by 



rup i^ik) = Re [Sp^h{k)] - k . 
For example, the mass of the particle mode obeys 



rrir, 



nir, 



rriT for — > 
2 rriT ioT k ^ T . 



(77) 



(78) 
(79) 



26 



The residue functions Zp{k,^) and Zh{k,^) govern the relative importance of particle 
and hole contributions to the thermal propagators. In the absence of interactions, one 
has Zp{k,n) = 1 and Zh{k,fi) = 0, in which case S'''{x,y; fj,) takes on the form given in 
Eqns. (llStiiyp . However, the departure from these limits is not perturbative in the strength 
of the interaction {e.g., the gauge coupling g), but rather depends strongly on the magnitude 
of the three-momentum, k. At = one has Zp = Zh = 1/2, while for k or order the 
temperature or thermal mass, one has Z^jZp ^1. In earlier work [1^, we studied the 
impact of the particle-hole structure of S^{x,y\\i) on fermionic contributions to the CP- 
violating source terms and CP-conserving relaxation rates generated by interactions with 
the spacetime varying Higgs vevs. We found that the gaugino and Higgsino contributions 
were dominated by large momenta (of order the gaugino and Higgsino masses), leading 
to negligible effects associated with the hole modes. In contrast, the hole modes generate 
non-negligible contributions to the quark relaxation rates. 

Here, we analyze the impact of hole modes on the three-body rates since it is not a 
■priori apparent that the loop integrals are dominated by momenta in the regime for which 
Zh/Zp <C 1. To that end, we consider one particular three-body process involving the 
massless fermion g, the massive fermion ^, and the scalar q. We compute the three-body 
rate vf"^ generated by the graph of FIG. [Ik but using Eq. (!73l) for the q Green's function. 
(Although we are choosing a particular interaction, this discussion is generic to any process 
involving one massless fermion, one massive fermion, and one massive scalar.) For illustrative 
purposes, we include only the gluino contribution, neglecting here those contributions from 
electroweak gauginos. 



20 F 




200 400 600 800 1000 
mq, (GeV) 

FIG. 2: We plot F^'^ (gluino-only) as a function of m^^. The solid and dotted lines denote, 

respectively, the particle and hole (xlO) contributions to this rate. The dashed line indicates F~'^ 
calculated using free thermal Green's functions with a thermal mass, corresponding to the last 
term in Eq. (j62b ). 

In FIG. 121 we show the resulting rrig-dependence of F^?'''^ for m-g = 200 GeV and T = 100 
GeV. The solid and dotted curves show, respectively, the particle and hole (xlO) contribu- 
tions to F~'^, calculated using the full quasiparticle Green's function given in Eqns. fl73H75l) . 
in the limit of zero thermal widths. The hole contribution is much smaller than the particle 
contribution because the process is dominated by momenta k ^ T; for these momenta, we 
always have Z^ < 10~^. 

For comparison, the dashed curve in FIG. m shows calculated using free Green's 

functions (|T9l) with the inclusion of thermal masses. The agreement between solid and 
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dashed curves in generally quite good; the majority of the discrepancy is due to the fact 
that Zp 1. At rriqj^ > 600 GeV, the rate T-.'^ is dominated by quasi-particles with 
momenta k > (few) x 100 GeV, for which Zp ^ 0.95; consequently, the agreement between 
solid and dashed curves is better than ~ 95%. For m^^ < 100 GeV, the rate is dominated 
by quasi-particles with momenta k ~ 100 GeV, for which Zp ~ 0.8; consequently, there 
is a 20% discrepancy between solid and dashed curved. In the present work, we compute 
r~''^ and all other three-body rates using free Green's functions with thermal masses, rather 
than the full quasi-particle Green's functions. Since these quasi-particle effects provide only 
small corrections to the three-body rate, we do not consider the more complicated case of 
two massless fermions and one massive scalar. 



3. TRANSPORT EQUATIONS AND Yb: ANALYTIC STUDY 

In this section, we present semi-analytical approximate solutions to the Boltzmann equa- 



we 



tions. In addition to the discussion that has already been presented in Refs. [23|, [24 
pay particular attention to the relaxation of particular charge densities towards superequi- 
librium in the diffusion wake ahead of the bubble wall. For this purpose, we first identify 
the condition that determines whether or not a particular interaction rate is sufficiently fast 
to lead to chemical equilibrium. We then show that fast supergauge interactions rates are 
a sufficient, but not necessary, condition for superequihbrium. Indeed, it turns out that the 
combined effect of Yukawa and triscalar interactions can also lead to superequihbrium for 
particular species. 



3.1. Conditions for superequilibrium 



Having specified the relevant interactions during EWB, we discuss here under which 
circumstances a particular interaction maintains chemical equilibrium and when this equi- 
librium is physically relevant. As an example, we consider supergauge interaction rates in 
the transport equations for q and q: 
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For the moment, we focus only on the gaugino interactions that maintain chemical equi- 
librium between q and q. Expressing these two equations in terms of chemical potentials, 
rather than charge number densities, and taking the difference, we have 
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Formally, chemical equilibrium corresponds to the equality of chemical potentials. 
Eq. fl82|) implies that the difference of chemical potentials, /ig — /ig, will relax to zero with a 
characteristic time scale 
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Consider now the presence of a density q {e.g. induced by decays W ^ q + t*) . In 
general, it will equilibrate with g on the time scale Tgq, so long as the corresponding reaction 
q + V ^ q is kinematically allowed. In the limit that gauginos become heavy compared to 
the temperature, we expect this equilibration process to be Boltzmann suppressed, that is: 

Teq — > CXD , niy OO . (84) 

This suppression of the equilibration rate arises as one would expect because the gauginos 
decouple from the plasma. On the other hand, taking the limit rrig oo leads to 

rif oc ^ . (85) 

This result is somewhat counter intuitive, since the g density becomes Boltzmann suppressed 
for large rrig. However, the corresponding chemical potential is proportional to q/kg, and 
kg also decreases with larger rriq. As a result, the chemical potentials of a particle and its 
superpartner may equilibrate quickly in the presence of gauginos, even when the sparticle 
mass is much larger than the temperature. This is because the sparticle chemical potential 
can adapt to the particle chemical potential by a small change in the physical sparticle 
density. 

From the standpoint of the generation of nieft, it is important to analyze the converse 
situation, namely, how the presence of a heavy superpartner affects the chemical potential, 
and hence the number density, of the lighter SM particles through the network of transport 
equations. From our numerical studies, we find that when a (s)particle [e.g. q), becomes 
heavy compared to the temperature, its transport equation effectively decouples from the 
system of transport equations. Consequently, even though a process of chemical equilibration 
involving this heavy (s)particle may take place quickly, its occurrence will be irrelevant for 
the generation of the lighter particle densities. As we will see below, this decoupling effect 
can lead to "bottlenecks" in a chain of reactions that might otherwise lead to chemical 
equilibrium and significant effects on particle number densities. 

With these considerations in mind, we will consider the equilibration time scales as in 
Eq. (1831) to be physically relevant only when the (s)particles involved are not too heavy 
compared to the temperature. In this regime, the corresponding fc-factors are never too 
suppressed with respect to their values in the massless limit. When more than two particles 
are involved in the equilibrating reaction, the corresponding equilibration rate will go like 

-e.^-rx5^1 . (86) 

For purposes of determining the criteria for various reactions to reach equilibrium on time 
scales short compared to other processes, we will take the longest possible time scale for 
each reaction by retaining only the minimum value of the 1 / kj appearing above: 

r^-i = r(^i'^2,...) I ^ I . (87) 



A sufficient condition for determining whether chemical equilibrium is maintained during 
the process of diffusion ahead of the bubble wall is then 



Teq < Tdiflf 



(88) 
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where 

Tdiff = 4 (89) 



is the diffusion time scale, with effective diffusion constant D defined below. If Teq > Tdiff, 
then an asymmetry between chemical potentials will diffuse ahead of the moving bubble 
wall faster than it will equilibrate away; therefore, equilibrium will be broken (unless there 
is another interaction fast enough to maintain equilibrium). When we evaluate all the 
interaction rates in Eqns. fl7Uati7Ur() numerically in Sec. HJ we will see which rates satisfy 
Tdiff /Teq > 1 SO that they maintain chemical equilibrium. (We note parenthetically that 
the Hubble rate if ~ 10 x T'^/Mpi is far too small to be relevant for EWB transport 



dynamics). In our previous papers [23|, |2J], we have taken = pC^i'^a,...) f^j. simplicity. 
The defintition flHTj) corresponds to a weaker requirement in the sense that Teq goes to zero 
in the case when all masses of the particles involved in the reaction are much larger than T, 
Now, one might expect that a sufficient condition for breaking superequilibrium is simply 
Eq. ([HI]). This expectation is false, as we can see from the following argument. Suppose all 
gaugino interaction rates vanished, while the Yukawa rates were infinitely large. Top quark 
Yukawa interactions in chemical equilibrium lead to the following relations among chemical 
potentials: 

(90a) 
(90b) 
(90c) 
(90d) 
(90e) 
(90f) 

Using these relations, one has 
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Note that the interaction leading to fl90bp is suppressed by the Higgs mixing parameter sin a. 



However, Eqs. fl9Tl) follow from Eqs. fl90|) even when Eq. fl90bp is excluded. Analogously, when 
the interactions mediated by bottom and tau Yukawa interactions are fast, it additionally 
follows that 

fJ'b = fJ'i , fi'e = fJ'i, fJ'T = IJ-T- (92) 

This implies that even in the limit rriy — > oo, where supergauge interactions are sup- 
pressed, superequilibrium can be maintained through Yukawa interactions alone. Fast su- 
pergauge interactions are a sufficient condition for superequilibrium, but not a necessary one 
in the presence of other interactions. Obviously, this argument does not apply to the first 
and second generation (s) quark and (s)lepton sectors, where the Yukawa rates are much too 
small to enforce superequilibrium in the absence of fast supergauge interactions. 

But also third generation (s) quarks and (s)leptons may not reach superequilibrium in the 
case when one or more of the (s) quark and (s)lepton species are very heavy. Then, particular 
Yukawa and triscalar rates can be small such that Teq defined according to Eq. (1871) fails 
the criterion Teq < Tdifr. As discussed above, this occurs because the heavy (s)quarks and 
(s)leptons can consitute a bottleneck for the transfer of the Higgsino charge to the non- 
supersymmetric particles. As a consequence, not all of the relations (l90l) need to be satisfied 
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at the same time and superequilibrium (19T]) does not follow. We encounter an example for 
such a situation in Section 14. 3[ 

To summarize, we identify two equilibration time-scales: 

• The time-scale t^^ in Eq. (183|) is useful to show that superequilibrium is maintained 
in the presence of light gauginos. In the case of a heavy q, the condition fig = fig is 
maintained due to the smallness of kg, but it is of little physical significance since the 
density of q is small. 

• The scale Tgq in Eq. (!87|) is useful to formulate sufficient conditions for equilibration. 
When Teq > Tdiff, there may or may not be a bottleneck preventing the establishment 
of chemical equilibrium on diffusion time-scales. Whether chemical equilibrium is 
maintained depends in this situation on the particular masses and interaction rates 
important for the network of reactions. 



3.2. Analytical approximation 

Before we explore numerical solutions to the diffusion equations, we briefly recapitulate 
some salient points of the analytical approximation that has been presented in 23|, l24 . 
First suppose that the relations fl91|92p hold, because of fast supergauge interactions or fast 
Yukawa and triscalar interactions. A relaxation of this assumption, among other things, is 
numerically studied in Section HI 

Then, it is convenient to introduce common particle and superparticle densities as 

X _x _ X _ _ 1,-1. 

— — — — — — , A — X + X , kx — kx -\- kx 1 

for X G {qi, Ui, di, i,r = r}, i E {1, 2, 3} ; (93a) 

H H H ~ 

=T7 = 7- , for i G {1, 2} , H = H + Hi + H2, ku = k^ + ku, + kn, ■ (93b) 
kui k^ kn 

Adding the equations for particles and their superpartners then immediately reduces the 
system of Boltzmann equations ( ffOl) to the more compact system of equations presented in 



Ref. [2J] that consists of six equations only. 

The next step towards an analytical solution is to realize that the sums of chemical 
potentials that multiply the Yukawa, triscalar^ and strong sphaleron rates can simultaneously 
be set to zero, as these rates are typically faster than the diffusion rate and the corresponding 
reactions reach equilibrium before densities can diffuse away. Explicitly, when imposing 

^ + ^-^ = ^-^-^ = ^-^-^ = 0, (94) 
kq kn kr kg ku ks ki kn kn 

in conjunction with approximate baryon number conservation (weak sphaleron transitions 
are out of equilibrium on diffusion time-scales r„s <^ '^difr); we can use these relations to 



^ We assume in this section that we are in a parametric region, where t-, b- and r- Yukawa and triscalar 
rates are fast. 
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eliminate all number densities except for H from the Boltzmann equations through 

Q = kqH = — H , (95a) 

kn kB + kQ + kr 

T = ktH = ~^ H (95b) 

^ kn ks + kg + kr 

B--sH^-'f .^l\^l\ H, (95c) 

kn ks + kQ + kr 

L = .,H = i,,tL—^s£^H. (95d) 

kn kiDi + kuUR 

Here, we have defined 

{ ; if << r: ■ 

If neither of these inequalities is amply fulfilled, the analytical approximation will not be 
accurate. Similarly, for ^ D^, Eq. fl95dp is in conflict with the fact that the densities 



R and L diffuse at a different rate ahead of the wall (unlike for the (s)quarks, which we 
assume here to diffuse at the common rate Dq, that is dominated by strong interactions). 
However, since ^ Dl, the the relation (]95dp is not a bad approximation compared to 



other uncertainties in the analytical calculation. This is because the right-handed (s)leptons 
diffuse on much larger distances away from the bubble wall, such that their local density can 
be neglected. Note that global lepton number conservation J dz{L + R) = holds, when 
neglecting weak sphaleron transitions. For more details on this point, see the discussion in 



Ref. 24 



Another feature of Eqs. fl95l) is that these relations imply that the contribution from third 
generation quarks to the factor that multiplies the strong sphaleron rate, Fss, vanishes. 
Since the CP-violating sources, triscalar couplings, and Yukawa interactions for the first 
two generation (s)fermions are highly suppressed by their Yukawa couplings, there exists 
no independent source for their densities apart from their coupling to the third generation 
via the strong sphalerons. Consequently, the vanishing third generation contribution to iVs 
implies that no chiral charge densities within the first two generations are generated 



This situation differs substantially from from that of earlier studies, where bottom Yukawa 



couplings are neglected [20|, |25l], and it leads to significantly different results for Yb and its 
dependence on the MSSM parameters. 

Applying the eliminations (p5|) . the resulting diffusion equation is 

v^H' - DH" = -VH + S , (97) 

and in the symmetric phase, ahead of the bubble wall, its relative accuracy is C(rdifr/ry). 
In this equation, the effective diffusion constant, source terms and damping rates are given 
by H 

D = Dn + DQ{-T-^B)+DL^Ln, 

r = — a T ) (98b) 

S = ^ ^ — — . (98c) 

1 + kt - Kb + VlHl 



32 



Using Eqs. (195|) The left handed fermionic charge density, that couples to the weak sphaleron 
is related to H as 

\kH ks + kg + kx ku kiDi + kuUR J 

In the symmetric phase, where the sources and relaxation terms are vanishing, the Higgs- 
Higgsino density is then given by 

H = ^e^™^/^ . (100) 

Assuming that the source terms are negligible for z < —L^/2, and that the relaxation terms 
have the particular form T{z) = 6{z)r, the normalization A can be found as 

dyS{y)^+ dyS{y)[-^ + |, (101) 

^ 



7± = i (v^ ± JvI + ATD] . (102) 



where 

= 2^ V ^» ' "~ J 

These analytic results, which have been discussed extensively in Refs. 23|, |2J|, suggest 
qualitative features one should expect from the full numerical solutions that we present in 
Section HI 

First, the presence of important bottom Yukawa interactions effectively quenches the 
contributions from the first and second generation quarks to nieft. This quenching arises in 
this regime because Yukawa-induced equilibrium involving third generation quarks leads to 
vanishing third generation chiral charge. Non-vanishing first and second generation quark 
densities arise only as required to maintain strong sphaleron equilibrium and thus, in this 
limit, also vanish. The resulting expression for nieft, given in Eq. (P^ ). depends only on 
quantities arising from third generation left-handed fermions. 

Second, the efficiency with which the non-vanishing H density - induced by the corre- 
sponding CP-violating source terms for Higgsinos - converts to nieft depends critically on 
the fc-factors associated with the right-handed third generation sfermions. The combination 
of Yukawa-induced equilibrium, superequilibrium, and local baryon number conservation 
implies that the third generation LH quark density induced by H depends on ks — krp. In 
the limit that the RH top and bottom squarks are degenerate, this contribution will vanish, 
while for a non-degenerate spectrum, the sign of this contribution will depend on which of 
the two squarks is heavier. The fc-factor dependence of the third generation LH leptons is 
more complicated since the LH and RH sleptons may diffuse at different rates. In the illus- 
trative limit of equal diffusion constants for the two chiral species, the lepton contribution 
depends on the geometric mean of the two fc-factors. If either the LH or RH sleptons become 
heavy compared to the temperature, the corresponding fc-factor is suppressed, signaling a 
decoupling of the slepton from the plasma and quenching the lepton contribution to nieft. 
Note that the sign of the lepton and bottom quark contributions are both opposite to that of 
the top quark contribution, implying that neglect of the former can lead to an overestimate 
of nieft ~ and, thus, of Yb - compared to the result when they are included. 

In the following discussion, we will see how these features emerge from the full numer- 
ical solutions to the coupled transport equations in regions of parameter space where the 
corresponding assumptions behind the analytic treatments are valid. We will also identify 
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TABLE I: Input parameters at the benchmark point. The masses for supersymmetric particles 
that do not occur in this table have been chosen to be 2TeV, such that they effectively deccouple. 

regions wherein these assumptions break down yet some of the quaUtative features persist, 
as well as regions wherein we find significant departures. 



4. TRANSPORT EQUATIONS AND Yb: NUMERICAL RESULTS 

We numerically solve the full system of Boltzmann Eqs. fITOl) in the presence of finite 
r^, Fy, Fss. Initially, we focus on one particular benchmark point, which is motivated from 
the requirement of a strong first order phase transition within the MSSM and for which our 
analytical approximations are valid (Section 14.11) . Then, we explore regions of smaller tan/3, 
where the analytical approximations break down because down-type Yukawa interactions 
do not equihbrate on diffusion time-scales (Section [4. 2p . We also investigate how robust the 
assumption of supergauge equilibrium is, in particular whether it can be maintained even 
when gaugino interactions are quenched (Section 14.31) . Technical details on how we obtain 
our numerical solutions are given in Appendix |X1 



4.1. Supergauge and Yukawa interactions in equilibrium 

First, we present a numerical solution for a particular point in parameter space. This 
point is chosen based on the following criteria: 

• it is motivated by existing studies of EWB within the MSSM, 

• there is a reasonable agreement between analytical approximations and numerical 
solutions, 

• the results share some qualitative key features with the scenarios recently presented 



in Refs. 23, 24 



The set of parameters that we choose is given in Table [B From these input parameters, 
we derive the parameters appearing in the Boltzmann equations f l70|) in the following way: 



The squark and slepton mass parameters in Table [I] are understood to be evaluated 
in the symmetric phase and without thermal corrections. We indicate this by the use 
of a capital M. The corresponding quark and lepton masses in the symmetric phase 
without thermal correction are zero. In addition, we derive the the masses of the Higgs 



34 



boson eigenstates from Mz and Ma as explained in Section [2l2| where we take account 
of the thermal corrections that are summarized in Ref. 24 



The source and damping terms are computed following 25, 2^, as explained in detail 
in Section 12.61 Additional input parameters here are the thermal widths of some 
particles X, which we denote as 7x. The particular values we adopt are again given 
in Table [H They are taken from Refs. 5^, 51 1 or motivated by the discussion therein. 



As mentioned above, we have chosen the mass parameters such that they are favorably 
disposed toward the viability of EWB in the MSSM. The light tji provides a strong first- 
order phase transition j3, 52, 5^. The heavy g is needed to increase the mass of the lightest 
Higgs boson beyond LEP II constraints. Making the first and second generation squarks 
heavy suppresses one-loop contributions to EDMs and precision electroweak observables. 

As a consequence of the large mass of g, there is wide agreement that EWB in the MSSM 
is viable only for CP- violation in the Higgsino/gaugino sector and only close the resonance 
region of Mi ~ (or M2 ~ and not from the quark/squark sector. Consequently, we 
choose Ml = \fi\ for our benchmark point to maximize the Higgsino CP- violating source. 
While we calculate our CP- violating sources as per Ref. |25|,|23, there is still some disagree- 
ment about the magnitude and parametric dependence of cf. [13, 31, 33|. However, 



in all of the present analysis, we keep fixed, as our emphasis is on the effects CP- 

conserving transport coefficients for a given CP-violating source. Therefore, our present 
work can be adapted to other calculations of simply by an appropriate rescaling, since 



all particle densities and the BAU scale linearly with 



In addition, following Refs. |30l . |5J] , we take the Higgs vev profiles to be 



v{z) = -vh{T) 



tanh I 2a 



f3{z) 



/?o(T)-^A/3(^l-tanh (s";^)) 



(103) 
(104) 



with a = 3/2, which provide an accurate analytic approximation to profiles obtained nu- 
merically in Ref. 13] • We follow Refs. (25I . 26 ] in our calculation of CP- violating sources 
and CP-conserving relaxation rates for quarks and Higgsinos. To facilitate comparison to 
other work, we provide the explicit numerical values for these results: 



H 

r(t,<?) 

^ M 



-9.4 GeV yism(l)^v^l3\z)v{zf 

1.4 X 10"^ GeV"^ X v{zf , 

6.2 X 10"^ GeV"^ x ylvu{zf , 

6.3 X 10"^ GeV"^ x ylv^izf , 
5.0 X 10"^ GeV"^ X y1vu{zf ■ 



(105a) 

(105b) 
(105c) 
(105d) 
(105e) 



As discussed earher, we neglect T^'lp and T^^^'^'^\ the relaxation rates for squarks and Higgs 



bosons, respectively; we defer a calculation of these rates to future study. Neglecting these 
relaxation rates will not have a large impact upon the final BAU to the extent that (a) 
our analytical arguments from Sec. [3] hold true, and (b) these rates are non-resonant and 
therefore much smaller than the relaxation rates for Higgsinos and quarks that we have 
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included. The rates F^^"^, F^^'^^-* affect the BAU only through F, which itself depends on 
the sum of all relaxation rates fl98bp : 



F oc F 



M 



M 



^ ^ H 



(106) 



Therefore, a proper inclusion of squark and Higgs boson relaxation rates would give only an 



H 



M 



correction to the analysis offered here. 

The degree to which we should expect an agreement between analytical and numerical 
solutions can be inferred from Table [Tll where we show each of the Yukawa and supergauge 
rates for our benchmark scenario. As we noted in Sec. [3l the appropriate sufficient condition 
that ensures that a particular rate leads to chemical equilibrium is idifr/req > 1. In the the 
second column of Table HTl we compute r^^is/Teq, where Teq is determined through dividing 



the corresponding rate from the first column by the appropriate fc-factor, as per Eq. 
The three body rates in Table [TTl are all calculated following the methods described in 



Sections 12.31 and 12.41 with one exception: F 



{t,q,H) 
Y 



After including thermal corrections, we 



have nit ~ ~ 65 GeV and niHi ~ 50 GeV, such that on-shell scatterings of these three 
particles are kinematically forbidden. Within the MSSM, we cannot alleviate this kinematic 
suppression, since a negative mass square for the Higgs boson needs to be present in the 
Lagrangian in order to ensure electroweak symmetry breaking. Since the rate Yy'"^'^^ plays 
a pivotal role in the network Boltzmann equations fl70|) that describe diffusion, and it is in 
general non-zero due to the off-shell effects and four body contributions, we treat it as an 
exception. From Ref. 



21 , we take 



<t,q,H) 



0.129^T, 

4:71 



(107) 



which is an estimate of the four-body contributions. We re-emphasize however that a more 
detailed analysis of the four body and off-shell contributions in the future would be desirable. 
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FIG. 3: Charge densities over z. Numerical results are represented by thick lines and analytical 
results by thin lines. Left panel: qs (pink, dot-dashed), £3 (green, dotted), nieft (red, solid). Right 
panel: H = Hi + H2 + H (blue, solid), qi + q2 (orange, dot-dashed). 



In FIG. 131 we display the numerical results for the number densities of selected species 
and compare them to the analytical predictions according to Section [31 Numerically, the 
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Rate 



V _ 

V _ 

V ^ 

V ^ 

V ^ 

p{5/) 

V 
p(T,r) 

-pit,Q,Hi) 
Z ^ 

p(6,Q,H2) 

y ^ 

X.~ 

■pit,Q,H2) 
Y 



r' 



Y 
Y 
Y 

-p{t,Q,Hi) 
^ Y 

-p{b,Q,H2) 
'- Y 

^ Y 

-pit,Q,H2) 
^ Y 

■p{b,Q,Hi) 

^ y 

Y ^ ^ 

Y 

Y 

^{T,i,H) 

Y 



Benchmark Value (GeV) 



5.3 
3.0 



3.3 
3.0 



1.6 
7.6 
9.3 



; X 10-1 
X 10^1 
1.8 X 10~^ 
; X 10^1 
X 10^2 

6.7 X 10-6 
X 10-1 
X lO-'^ 
X 10-^ 

1.8 X 10-s 



3.3 
4.2 



X lO-'^ 
X 10-^ 



8.1 X 10-^ 

8.0 
5.9 X 10-1 
8.9 X 10-2 

2.3 

9.6 X 10-1 
1.4 X 10-1 
1.8 X 10-2 



3.8 X 10-^ 

3.3 X 10-5 

4.4 X 10-6 

7.9 X lO-'^ 

3.7 X 10-1 



Tdifi/T( 



cq 



3.4 X 10 
2.6 X 10 

6.4 X IQ-^^ 
1.6 X 10 

2.1 
6.7 X 10-4 

2.8 X 10 
3.6 X 10-5 

2.1 X 10-5 
3.0 X 10-6 

1.6 X 10-5 

2.7 X 10-6 

5.2 X lO-'^ 

2.8 X 10^ 

2.1 X 10 
7.6 

8.2 X 10 
3.4 X 10 
1.4 X 10 

6.5 X IQ-i 



2.4 X 10-3 

2.3 X 10-3 
3.0 X lO-'' 

10-5 
10^ 



,8 X 



6. 
1.3 



X 



TABLE II: Yukawa, strong sphaleron, and supergauge interaction rates for the benchmark param- 
eters and comparison to the diffusion time scale. 



diffusion time scale is 



D 



TdifT 



2.0 X 10^ GeV" 



10-2^ s. 



(108) 



Since Tdis/Teq > 1 for those rates in Table HTl that do not involve the ultraheavy sfermions 
g and £ and that are not suppressed by sin a, our numerical solutions match our analytical 
expectations to a good extent. Regarding the analytical solutions, the following additional 
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comments are in order: 

• The overall normalization of the analytic estimate is slightly larger than for the nu- 
merical result. This is because the damping rates in the broken phase are larger than 
the equilibration rates, and therefore the error of the analytic approximation is not 
under control for z > 0. 

• Sufficiently far ahead of the bubble wall, the prediction for the ratios of the particular 
charge densities is accurate. In this region, the analytical approximation is well under 
control. However, nieft close to the bubble wall contributes significantly to Yg, and 
in this region the analytic approximation is less reliable. Close to the bubble wall 
and within the bubble, the relevant time-scale is no longer Tdis in Eq- fl89l) . but F"^ 
[cf. Eqs fl98bt llOll 11021) ]. Comparison of the numerical values for the relaxation 
rates (11051) to the equilibration rates in Table [ITl shows that Ty,TY -C TMy^H, such 
that the analytical approximation is not justified in that regime. In particular, since 
the densities q and i have opposite sign [as is generic for > m^, cf. Eq. fl^ ]. the 
analytic result for riien is rather inaccurate. Here and in the following, we therefore 
refrain from a direct comparison of the analytical predictions for the baryon asymmetry 
to the numerical answer. However, the analytical formulae are still very useful for 
understanding the behavior of the numerical solutions qualitatively. 

• The results in FIG. [3] do, indeed, refiect many of these qualitative features. In par- 
ticular, in the right panel, we observe that the total first and second generation LH 
quark + squark densities are negligible compared to the third generation LH quark 
+ squark and RH tau + stau densities shown in the left panel. This feature follows 
from the approximate vanishing of the third generation contribution to N^, leading 
to the near absence of any induced first and second generation densities, as explained 
above. In addition, the relative signs of the H, and £3 densities ahead of the wall 
follow closely the expectations based on Eqs. (p5!) (recall that < kx for our fiducial 
parameter choice). 

Besides, due to the choice of parameters in Table [H the scenario presented in this Section 



shares the key features with those that are presented in Refs. [23|, |2J]. In order to discuss 
these features, it is instructive to consider certain combinations of chemical potentials, that 
are displayed in FIGs. IH E] under the label "tan /5 = 15" . The following points are of 
importance: 

• Superequilibrium is maintained on diffusion time-scales Fjj^, as it is exhibited by the 
fact that the chemical potentials for particles and their superpartners are identical 
sufficiently far away from the bubble wall, cf. FIG. H] ("tan/3 = 15"). This feature is 
crucial for the validity of the analytical approximation in Section 13. 2[ The assump- 
tion of superequilibrium has also been used in deriving the reduced set of Boltzmann 
equations, which is presented in Ref. [2^ and follows from Eqs. fl70|) . 

• In deriving the analytical approximation, it is assumed that the relaxation rates Fy 
are fast compared to the diffusion time-scale Fjifj. This implies that the combinations 
of chemical potentials that multiply these relaxation rates should vanish. From FIG. E] 
("tan/3 = 15"), we see explicitly that Higgs bosons and Standard Model fermions 
satisfy the corresponding equilibrium conditions. We have checked that the same is 
also true for the additional supersymmetric particles. 
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The fact that Yukawa and triscalar interactions involving (s)bottoms and (s)taus are 
in equihbrium leads to an important change in the flavor dynamics ahead of the bubble 



wall, that has not been appreciated in the literature |20l . |25| before our resent work [23|, 
13]. First, since we have a moderately large value of tan/3, also lepton and slepton 
densities are induced and equilibrate ahead of the bubble wall and give an important 
contribution to electroweak baryogenesis. The same is true for the (s)bottom-particles, 
as they have larger Yukawa and triscalar interactions than the (s)tau. Due to the 
presence of the strong sphaleron, the equilibrium of bottom Yukawa and triscalar 
interactions has an additional consequence: The combination of chemical potentials 
A^5 that multiplies in the Boltzmann equations flTOj) is vanishing for zero densities 
of flrst-generation quarks. Indeed, from FIG. [3] we see that no asymmetry in flrst 
generation quarks is diffusing ahead of the bubble wall. 



4.2. Dependence on tan /? 

We now take the same parameters as given in Table [T] for our example point, but consider 
values for tan (3 G [1.5; 20]. The effect of this on the resulting baryon asymmetry is displayed 
in FIG. [6|, where we display the ratio of Yb obtained from the numerical simulations and 
the observational value Y^^^^ . We see that Yb/Y^^^^ increases for smaller values of 
tan f3. This is because smaller values of tan f3 imply smaller down-type Yukawa couplings. 
Therefore, a smaller lepton density is generated ahead of the wall, as can be seen when 
comparing FIG. [7] to FIG. [31 Since quark and lepton asymmetries contribute with opposite 
sign (provided m^; > mj, as it is the case here), small values of tan/3 lead to a weaker 
cancellation in the left-handed fermion density and therefore a larger baryon asymmetry. 

Besides, FIG. O ("tan/? = 1.5") exhibits that for tan/? = 1.5, neither the interactions 
mediated by yb nor t/r maintain equilibrium. However, even for values of tan/3 close the 
lower bound that is theoretically allowed, a non-negligible density in 6-quarks is produced. 
Apparently, the analytical approximation is only reliable for tan/? > 15, when all third 
generation Yukawa interactions are in equilibrium on diffusion time-scales. For tan /3 < 15, 
the interactions mediated by i/r are out of equilibrium on diffusion time-scales, and for 
tan/3 < 5 the same is true for the interactions mediated by yi,- Yet, non-negligible densities 
of down-type fermions can occur in general ahead of the bubble wall. We also note that 
since for tan/3 = 1.5, bottom (s)quarks do not equilibrate ahead of the wall, a non-negligible 
density of flrst-generation quarks is generated, cf. FIG. [71 This density of left-handed quarks 
of the flrst two generations is opposite to the third generation density, and it has therefore 
the effect of suppressing the baryon asymmetry {cf. the graph without taking account of 
leptons in FIG. [6|). For large yb, the equilibrium of axial charges that is maintained by the 



strong sphaleron is satisfled by the relation fib + IJ-t — fJ-q = and fig^ ^ = |23|, [2J], whereas for 



negligible yb, one flnds Q12 = 2(Q + T) (where Q and T are of opposite sign and |T| > \Q\), 



hence and Hg^ ^ being opposite [20|, |25|. Note that the analytic formulae presented in 



Refs. (20I . 25| are not applicable for tan/3 = 1.5, since even though bottom quarks do not 
completlely equilibrate, a sizeable density of them is yet present ahead of the wall. 

Comparison of the graphs with and without leptonic densities taken into account in 
FIG. [6l we also observe that for values as small as tan/3 = 1.5, there is still a sizable 
leptonic contribution to nieft. We also see that fast interactions mediated by Binos and 
Winos still maintain superequilibrium ahead of the bubble wall, as exhibited in by the small 
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FIG. 4: Chemical potentials over z illustrating supergauge-equilibrium. Key: (pink, dotted), 
p,H2 (green, dashed), /ijj (red, solid), Hq,i,t,b,T (red, solid), fJ'gjj;h,T (P™1^' dotted). 
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FIG. 5: Chemical potentials over z, illustrating Yukawa-equilibrium. The key is fit — /^g, fJ-q — Mfe) 
fi£ — fir (red, solid), fiHi, f^H2 (green, dashed). 
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FIG. 7: Charge densities over z, at the fiducial point but with tan/3 = 1.5. Left panel: qs (pink, 
dot-dashed), £3 (green, dotted), nieft (red, solid). Right panel: H = Hi + H2 + H (blue, solid), 
Qi + Q2 (orange, dot-dashed). 

tan/3 region shown in FIG. H] ("tan/3 = 1.5"). 

We now investigate the impact of the finite rate of supergauge interactions. For that 
purpose, in FIG. [6l also the results of a simulation where superequilibrium is enforced 
(leading to /i^ = everywhere) are displayed. The corresponding plots of the chemical 
potentials can be seen in FIG. [8] (Fy 00). At large tan /?, compared to the case with finite 
supergauge interactions [FIG. S] (tan /3 = 15)], /i^ is enhanced, since it adapts to some extent 
to iJ,j [cf. FIG. m (tan/3 = 15)]. On the other hand, since g is superheavy {i.e. the density 
of q is very small), there is no corresponding enhancement of /ig. In addition, imposing 
superequilibrium slightly favors the production of leptons close to the bubble wall due to 
the summation over the various supersymmetric production channels, whereas the quark 
production rate is already comparably large. The combination of these effects leads to an 
enhancement of the chiral lepton asymmetry which in turn suppresses Yb- For small tan/3 
the situation is more complicated, since now also first generation quarks play a role. Since 
for finite supergauge interaction rates superequilibrium is still violated to some extent close 
to the bubble wall, the apparent agreement of the curves with finite and infinite supergauge 
interaction rates in FIG.lHlfor small tan/3 is accidental. Note also that for large tan/3, since 
the contributions of i and q to njeft have opposite sign, the inaccuracy incurred by assuming 
finite supergauge interactions is substantial, up to the extent that the predicted value of Yb 
can flip sign. 



4.3. Absence of supergauge interactions 

We now investigate in more detail how superequilibrium can be maintained even in ab- 
sence of supergauge interactions, as discussed in Section 13. 1[ We first note that since we 
have taken its mass to be 2 TeV, the interactions of the gluino are suppressed to an extent 
that they are negligible. Now, we set in addition all interactions F^ appearing in the Boltz- 
mann equations (1701) to zero. While this procedure could also be mimicked by taking the 
Bino and Wino masses to be very heavy, we note that for electroweak baryogenesis within 
the MSSM, it is at least required that either Mi ~ /i or M2 — /i, in order to have resonant 
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FIG. 8: Chemical potentials over z illustrating supergauge-(non-)equilibrium for the case when 
supergauge interactions are infinitely slow and infinitely fast, respectively. Key: (pink, dotted), 
IJ.H2 (green, dashed), fi^ (red, solid), fiq,e,t,b,T (red, solid), fJ'qjj^^:^: (pink, dotted). 



CP-violation and to produce a large enough baryon asymmetry. Therefore, the limit taken 
in this section may be considered as a theoretical exercise. On the other hand, it may be 
conceivable that M2 ^ 1 TeV and Mi ~ fi, but also fi ~ Mi ~ ^tbr- such a case, 
three-body interactions between the Higgsino and right handed fermions are kinematically 
not allowed at zero temperature. Besides, if the CP-violating source in supersymmetric 
scenarios different from the MSSM is not originating from Higgsino-gaugino mixing, there 
may be no obstacle for successful EWB with heavy gauginos. 

For the purposes of this example, we again take the parameters from Table [U but we set 
Mj = 100 GeV. For which species superequilibrium is maintained and for which it is broken 
can now be inferred from FIG. H] (F^ = 0; Mj = 100 GeV): while {q, q} and {t, t} do not 
satisfy superequihbrium, {i, £} and {r, r} do. 

To give an explanation of these observations, we first note that since q is superheavy, 
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FIG. 9: Charge densities over z, at the fiducial point but without supergauge interactions. Left 
panel: qs (pink, dot-dashed), £3 (green, dotted), nieft (red, sohd). Right panel: H = Hi + H2 + H 
(blue, solid), qi + q2 (orange, dot-dashed). 



the chain of equilibrium conditions (pOj) is broken. To see this, we note that fiq is sizeable, 
even far ahead of the bubble wall. However, the number density q is small, since we have 
taken the left-handed squark to be superheavy, rriq = 2 TeV. This implies that for example 

Fy^'^^/k^ <^ Fdifj and Vif'^^ /kn^^ <^ Tdifr- Therefore, a sizable value of /Xg does not need 
to enforce a large density of Hi , for the simple reason that the physical density q is small. 
In FIG. [5] {Ty = 0; M-jr = 100 GeV) it is exhibited, that even though superequilibrium is 
violated, Yukawa equilibrium is still intact for Standard Model fermions and Higgs bosons 
on diffusion time-scales. _ _ 

In contrast, in the down- type sector, b, r and i are not heavy compared to T, such 
that these particles can mediate the equilibration of H2. Consequently, superequilibrium is 
maintained here on diffusion time-scales according to the argument given in Section 13. 1[ 

In FIG. |6] (green dashed curve), we also show a simulation for the parameter set as in 
Table [I (but now again with Mj= 2 TeV), but with all supergauge interaction rates set to 
zero. The profiles of chemical potentials are displayed in FIG. |8] (F^ 0; Mj = 2TeV). We 
find that Yb is enhanced when compared to the cases with finite or infinite F^. Comparing 
to the plots in FIG. [8] {Ty —> 00), we see that in the absence of supergauge interactions due 

to the q and £ bottlenecks, superequilibrium is broken in both the up and the down type 
sector. In addition, we observe: 

• The signs of the t and r chemical potential are now reversed. 

• The magnitudes of the t, b, and £ densities are reduced. 

• The magnitudes of the Hi^2 are significantly suppressed. 

• The magntiude of the q density increases. 

On general grounds, the absence of superequilibrium implies a degrading of the overall 
efficiency with which H density (induced by the CP-violating source) is transferred into the 
SM fermion densities. At the same time, the detailed balance between these densities and 
their net effect on Yb change substantially. Without a robust analytic framework for treating 
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this case, we can only speculate on the reasons why these changes result in an increase in 
Yb- Nonetheless, it is clear that the complete decoupling of supergauge interactions from the 
transport dynamics can have a substantial impact on the predicted baryon asymmetry. This 
siutation may be particularly relevant to extensions of the MSSM that can accommodate 
sizable CP-violating sources and heavy gauginos. 



5. CONCLUSIONS 

In this work, we have generalized existing approaches to EWB to account for finite 
supergauge interaction rates. We have developed numerical solutions to the resulting Boltz- 
mann equations that describe the diffusion processes. For particular illustrative points in 
parameter space, we have presented numerical solutions. When superequilibrium holds 
and all Yukawa interactions fully equilibrate, these examples agree with results published 



earlier [23|, |2J]. In turn, in the absence of superequilibrium or when down-type Yukawa 
interactions only partially equilibrate, our solutions show sizable deviations from earlier 
results. 

Regarding the consequences of finite supergauge interactions for supersymmetric EWB, 
our conclusions follows: 

• In models for EWB with light gauginos [e.g. with masses not much heavier than T), 
such as the MSSM with a Higgsino-gaugino CP violating source, superequilibrium is 
a robust assumption ahead of the bubble wall. 

• In models with heavy gauginos - such as extensions of the MSSM that do not require 
light gauginos for the CP- violating sources - superequilibrium may be restored through 
the network of Yukawa and triscalar interactions. However, this chain of reactions may 
be broken when one of the superpartners becomes heavy compared to the temperature, 
thereby leading to a bottleneck. We have presented an illustrative example of this 
situation, where in particular the up- (s) quark sector violates superequilibrium, whereas 
superequilibrium is maintained for down-type (s)quark and (s)leptons, cf. FIG. H] 
(F^ = 0; Mj= 100 GeV). 

• More generally, the assumption of exact superequilibrium or the neglect of supergauge 
interactions may lead to significant errors in the prediction of Yb, as illustrated, re- 
spectively, by the yellow and green curves in FIG. HI In the former case, which may 
apply in models with heavy gauginos, the assumption of superequilibrium may lead 
to larger left-handed particle densities than actually occur as one forces them by hand 
to match the corresponding sparticle densities. As indicated by the yellow curve of 
FIG. m the effect of this error is amplified because the LH squark densities are small, 
the LH slepton densities are relatively large, and the corresponding particle densities 
contribute to nieft with opposite sign. As a result, the enhanced negative lepton con- 
tribution suppresses the baryon asymmetry. The green curve illustrates the converse 
dynamics, wherein the neglect of supergauge interactions leads to a lepton contribution 
that is smaller in magnitude with a correspondingly larger Yb- 

• When mj < (as it is typically the case in the MSSM with a strong first order 
phase transition), the leptonic contribution to nieft is opposite in sign to the one from 
quarks. The leptonic contribution is suppressed for small values of tan/?, but yet needs 
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to be taken into account for values as small as tan (3 = 1.5 (c/. FIG. |6]). In order to 
obtain a large asymmetry in the MSSM, small values of tan/5 are therefore needed 
from the solution of diffusion equations. This points into the same direction (but it 
is a different effect) as the suppression of the CP- violating source for large values of 



tan/3 [37|, 153|. Note that we have not taken the impact of tan/5 on the CP- violating 
source into account in this work, in order to disentangle this effect from the diffusion 
solutions. However, large values of tan (3 are yet interesting for EWB in the (M)SSM, 
since the Bino phase can be of order one from present EDM limits, which allows for a 



value of Yb in that region of paramater space |42|, |43 



From our results, we can draw the following conclusions regarding the relying on nu- 
merical solutions to the diffusion equations when one is interested in quantitatively reliable 
predictions for Yb'. 

• The numerical solutions are accurate close to the bubble wall and within the bubble, 
where the analytic description is not under control. This may have a sizable impact 
on the result for the BAU, which is obtained from integrating over nieft. Even in the 
presence of gauginos that are not heavy compared to the temperature, it is advisable 
not to impose superequilibrium, as it is usually broken close to the bubble wall. This 
is of particular importance when particular densities that contribute to nieft cancel, as 
it is often the case. 

• For tan/? < 15, r- Yukawa couplings or both bottom- and r-couplings equilibrate in- 
completely, but there is always a non-negligible density of b and r. Since there is no 
analytical method yet in order to describe this partial equilibration, the numerical 
solution is necessary for predictions in this parametric region, which is of particular 
interest for EWB in the MSSM. 



In closing, we emphasize that our conclusions apply to any scenario of supersymmetric 
EWB, beyond the MSSM. Our discussion has been nearly independent of the exact nature 
of the CP-violating source, which may arise from Higgsinos, squarks, or an extended sector 
not present in the MSSM. What is certain, however, is that in order for EWB to work, 
CP- violation must be communicated to the left-handed matter fermion sector, which in 
general will give rise to CP- violating asymmetries for squarks, quarks, sleptons and leptons 
of various flavors through supergauge, Yukawa, triscalar and strong sphaleron processes. 
The complete spectrum of gauginos, squarks and sleptons as well as the value of tan P are 
relevant for the precise determination of the baryon asymmetry in all supersymmetric EWB 
scenarios. 
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APPENDIX A: NUMERICAL METHODS 

We have developed two independent numerical codes to solve the system of Boltzmann 
equations. The first method, a modified version of the Euler method, is best illustrated for 
the differential equation 

y"{z) + a{z)y\z)+b{z)y{z) = s{z) , (Al) 

where a{z), b{z) ~ constant and s{z) ~ far enough away from z = 0, say, for z < zi < 
and z > Z2 > 0. We exploit the fact that Eqn. ( 1A1[) is homogenous and exactly solvable for 
z < Zi and z > Z2. Therefore, we can write 

, , f Ae^-"" z < zi , 
[ Be + z > Z2 

where X± > 0. Because a, b are known functions, X± are known as well; however. A, B 
are unknown coefficients. The characteristic equation resulting from Eqn. ( lAll) permits 
additional roots, which lead to solutions with diverge at z ^ ±oo; we implement the 
boundary conditions that 

lim y{z) (A3) 

by discarding these solutions. The next step is to use the Euler method to determine y{z) 
in the interpolating region where zi < z < Z2. We discretize this region into steps with 
length A. We begin with 



y{z,) = Ae'~^^ (A4) 

y'{zi) = AA_.e^-^i (A5) 

y"iz,) = AX'_e'-^^ ; (AG) 

then we iterate forward: 

yizi + A) = y{zi) + y\z^) A + ... (A7) 

y'(^i + A) = y'(^i)+/(^i)A+ - (A8) 
y'\zx + A) = -a(zi + A) y\z^ + A) - 6(2:1 + A) y{z^ + A) + s{zx + A) , (A9) 

where the denotes the possibility of including higher order terms if needed. Ultimately, 
after iterating from z\ to z\ + A^A = ^2, we obtain 

y{z2) = fi{\-,z,)A (AlO) 

y'{z2) = f2{X-,Z,)A, (All) 
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where /i^2 are simply numbers (which depend on A_ and zi). Finally, we merely have to 
solve the equations 



-A+22 



yiz2) = M\^,z,)A = Be 
y\z2) = f2{X-,z,)A = BX+e 



-X+Z2 



(A12) 
(A13) 



to determine the unknown coefficients A and B. The application to the system of coupled 
densities fl70atl70rl ) follows by generalizing Eqn. (lAll) to a matrix equation for the vector of 
densities y = (t/j, t^, ...) and performing some diagonization gymnastics. 

The second numerical way we employ to solve the diffusion equations is a relaxation 
method, see e.g. [55|. We decompose the diffusion equations into 2A^ ffist order differential 
equations 

yl{z) + ji,{z)yj{z) = Si{z) , (A14) 

where '~fij{z) depends on the interaction rates, the wall velocity and the diffusion constants 
and Si{z) on the source. Among the yi, N components represent the charge densities and 
A^ the derivatives of these densities with respect to z. 
Then, we discretize (1A14I) for M interior points as 



z'^ + z'' ^ 



y^ + yr 



z'' + z'' ^ 



0, (A15) 



where Zk = ^min + (^max — Zram)k/M. Two additional sets of equations at the exterior points 
Zmin ^ and ^max ^ foUow from the boundary conditions. Here, we impose that the 
charge densities are vanishing far away from the wall, that means for some pair of large 
negative and positive values of 2;min and -Zmax- 

We start with a initial guess yf = for all k and i. The yf = are then updated by 
solving the linearized approximation to (lAlSp 



E^{y'' + Ay',y'- 



Ay 



fc-l^ 



EHy^y 



k-l\ 



+ 



dEf{y^,y^- 



-Ay^ + 



dy 



k-l 



dE'l{y^,y^) 
dyj 



(A16) 

subsequently for Ay^~^ and Ay^. An improved approximation to the solution is then given 
by y'^ ^ + a Ay'', where a is a positive constant of order one, to be chosen such that fast 
convergence is achieved. For the present problem, a relative accuracy of one part in 10^° is 
typically attained after two or three iterations. 
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